Gap formation and soft phonon mode in the Holstein model 
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We investigate electron-phonon coupling in many-electron systems using dynamical mean-field 
theory in combination with the numerical renormalization group. This non-perturbative method 
reveals significant precursor effects to the gap formation at intermediate coupling strengths. The 
emergence of a soft phonon mode and very strong lattice fluctuations can be understood in terms of 
Kondo-like physics due to the development of a double-well structure in the effective potential for 
the ions. 
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Despite the many years of study of the electron-phonon 
interaction in metallic systems, there remain fundamen- 
tal problems that have yet to be resolved; particularly 
in the strong-coupling regime and in conjunction with 
strong electron-electron interactions. A solution to these 
problems will be required to understand fully phenomena 
such as the colossal magnetoresistance effect in mangan- 
ites M. Also in the metallic alkaline- doped Ceo-based 
compounds, high critical superconducting temperatures 
have been observed 0. These materials are known to 
have strong electron-phonon and electron-electron inter- 
actions Ejj. Recent photoemission experiments indicate 
strong electron-phonon coupling in the cuprate high- 
temperature superconductors ||. There is a clear need 
of theoretical techniques to tackle these problems in the 
strong-coupling regime. 

Although electron-phonon problems involving one or 
few electrons can be solved to very high accuracy |Q, |J, 
so far there are no comparably accurate approaches for 
the many-electron case relevant to metallic systems. In 
this letter we study the simplest realization of electron- 
phonon coupling: The Holstein model with finite electron 
density describes the coupling of Einstein (LO) phonons 
to the density of electrons of a non-degenerate conduction 
band: 
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No general exact solution of this model is known for sys- 
tems with finite electron density, even in the limit of 
infinite spatial dimensions (d = oo). This limit takes 
local quantum fluctuations fully into account and has 
proved to be a powerful tool in understanding strongly 
correlated systems |], fjj. Although exactly solvable for 
d = oo, the case of a single electron in the band || 
is physically very different from the many-electron case 
since no electron-electron pairing (bipolaron formation, 
superconductivity etc.) can occur. Another, more in- 
structive limiting case is the static limit coq = 0, where 



the phonons are replaced by a static displacement of the 
lattice ('static' or 'classical' approximation). Extensive 
calculations in this limit for d = oo have been presented 
in Ref. B, However, it is immediately clear that this 
static limit neglects all possible effects stemming from 
the quantum nature of the lattice excitations. In the op- 
posite limit of coq — > oo the lattice reacts instantaneously 
to the state of the electrons. Here, the Holstein model 
can be mapped onto a non-retarded attractive Hubbard 
model H by integrating out the phonons. The Hub- 
bard model has been intensively studied, and much re- 
cent progress has been based on using the d = oo limit. 
Although the large-wo limit is not physically relevant, it 
is still a useful point of reference for getting an overall 
understanding of the physics of the model. Of physical 
concern for applications are relatively small phonon fre- 
quencies of the order of ujq ~ 0.01W — 0.2W (W is the 
width of the electron band). In this regime, the Migdal- 
Eliashberg diagrammatic approach has been used p0| . 
The main feature of this approach is the neglect of ver- 
tex corrections. A sufficient criterion for its application 
is usually ujq/W <C 1. However, at least for half-filling 
there is evidence that this approach breaks down for in- 
termediate coupling strengths g |l(], |ll| . There have also 
been a number of perturbative schemes going beyond the 

^ Migdal-Eliashberg approach and including some vertex 

Corrections §, 0,P,0. 

In the dynamical mean-field theory (DMFT) a lattice- 
model is mapped onto an associated impurity model. The 
parameters of the associated impurity model are related 
to the Green's function of the lattice model by a self- 
consistency condition. This mapping becomes exact in 
the limit of infinite spatial dimensions (d — oo). The 
method is described in detail in Ref. |7|. One of the most 
precise techniques for solving the associated impurity 
model for low temperatures is the numerical renormal- 
ization group (NRG) H|, [L7j . It is capable of resolving 
very low energy scales, and gives information about the 
excitation spectrum over the whole real-energy axis. The 
combination of DMFT and NRG has helped to solve a 
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FIG. 1: Electronic spectral density p{uS) = — ilm((cio-; cJ CT )) 
from the DMFT-NRG calculation for ujq = 0.05 and various 
coupling strengths g. 



number of open questions regarding the Mott transition 
in the Hubbard model J?], [HI . The application of the 
DMFT to the Holstein model ([[]) leads to the Anderson- 
Holstein impurity model, which is essentially an Ander- 
son model with additional coupling of a local phonon 
mode to the impurity site: 
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An extensive study and discussion of this model is pre- 
sented in Ref. ^0|, which also gives details of the general- 
ization of the NRG to this situation pj^ . In this letter, 
we present and discuss results obtained for the (lattice) 
Holstein model (|l|) using the NRG in conjunction with 
the DMFT. 

We calculate the single-electron spectral function self- 
consistently within the DMFT approach. The local 
phonon propagator d{uj) = ((bi,b\)) for the Holstein 
model as well as spin- and charge susceptibilities for the 
associated impurity model can easily be obtained by this 
method. The phonon propagator d(w) can also be cal- 
culated from the charge susceptibility of the impurity 
model: 

d(u) = doM+ ff 2 doM 2 X cM, (3) 

where do(uj) = (u> — luq + i0 + )~ 1 is the phonon propaga- 
tor for .g = 0, and Xc( w ) is the charge susceptibility of 
the associated impurity model. Another frequently dis- 
cussed phonon propagator, D(lu) = (({hi + b\), (bi + b\))} 
can be calculated by a similar formula, which is obtained 
by replacing do(uj) by Dq(cu) = 2uiq/(uj 2 — ui 2 , + i0 + ) in 
Eq. (|. 
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FIG. 2: Quasiparticle weight Z = (1 asrM"=o) 1 as func- 
tion of g from the NRG and the Migdal-Eliashberg (ME) cal- 
culation 



In the following numerical results, we chose for the 
uncorrelated system [g = 0) a semielliptic density of 
states for the conduction band. Its width W = 1 de- 
fines the energy unit used throughout this paper. We 
only consider the particle-hole symmetric case and, in 
analogy to the Mott-transition in the Hubbard model in 
d = oo, suppress long-range order, which corresponds 
to suppressing anti-ferromagnetic order in the Hubbard 
model. Unless otherwise noted, the phonon frequency is 
taken to be luq = 0.05. All calculations were performed 
for T = 0, but our method can be extended to finite 
temperatures [ |l9| . 

In Fig. |l|, the electronic spectral function is plotted for 
various values of g with uj = 0.05. For weak coupling, a 
small feature emerges at the Fermi energy (ui = 0). With 
increasing g, this peak becomes narrower and more pro- 
nounced. This behaviour is qualitatively similar to that 
found within the Migdal-Eliashberg (ME) approach , 
the quantitative difference is the enhanced narrowing in 
the NRG calculations. 

At intermediate coupling, the central feature becomes 
very narrow, and two broad peaks emerge above and be- 
low the Fermi energy. These are entirely absent in the 
ME approach. At some critical coupling g c w 0.099, the 
central peak disappears and a gap opens between the two 
upper and lower bands [^2| . For all g < g c , the system 
is a Fermi liquid with ImE(tj) ~ to 2 for small energies. 
This manifests itself in Fig. |lj in the pinning of the spec- 
tral function at the Fermi energy. 

The quasiparticle weight Z = (1 — ^§^U=o) _1 is 
shown in Fig. || as obtained within the NRG and the 
ME calculation. In both cases, Z decreases with increas- 
ing g. Up to g « 0.05 both lines coincide, but then the 
NRG curve decreases faster. The ME calculation breaks 
down before Z reaches 0.1. 

We take a closer look at the large ojq limit in Fig. ||, 
where the electronic spectral function for l>q = 1 and 
3 are plotted for three values of g well below, close to, 
and above g c . The two broad bands discussed above for 
luq = 0.05 split now into two each. For luq — 3, the 
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FIG. 3: Electronic spectral function as in Fig. |y, but for larger 
ujo = 1 and 3. The values of g are given in units of g c which 
is g c = 0.69 (1.48) for lo = 1 (3). 



higher-energy ones are not visible any more on the plot- 
ted scale. These bands are multiphonon bands and are 
shifted approximately by ioq. The lower-energy ones are 
bipolaron bands located at the bipolaron binding energy 
A = 2g 2 /luq. The weight of the multiphonon bands van- 
ishes with increasing loq ]2(]] , and for sufficiently large loq , 
they can be neglected. The remaining excitation spec- 
trum corresponds to that of an attractive Hubbard band 
where the effective interaction \U\ — A corresponds to 
the bipolaron binding energy. 

The Mott metal-insulator transition in the Hubbard 
model shows a parameter regime with U c ± < U < U c = 
U C 2 where metallic and insulating solutions co-exist ('hys- 
teresis'). For large ojq, the effective |E4i,2| = 2g 2 12 /too of 
the Holstein model coincide with the values known from 
the Hubbard model |lS|| . For smaller u>o, the hysteresis 
region shrinks, and finally, for uj$ = 0.05 no hysteresis is 
detectable (g c i — g c2 ). 

To gain further insight into the results presented so 
far, let us look at the mean- field solution of model ([!]). 
The classical field x is self-consistently determined as 

x = vkS b + 6t > = ~ §)>• For sma11 9 

the excitation spectra remain unchanged from the .9 = 
case since any effects of the distortion are cancelled out by 
a change in the chemical potential. The system becomes 
unstable towards charge-order at a critical coupling gi mf 
which for luq = 0.05 is <?e mf ' = 0.085. If one (artificially) 
restores the symmetry, one obtains an electronic excita- 
tion spectrum consisting of 2 peaks shifted by gx m { above 
and below the Fermi energy. The self-consistently calcu- 




FIG. 4: Lattice fluctuations {x 2 - (x) 2 ) for NRG and Migdal- 
Eliashberg (ME) calculation. Additionally the value of x^f 
as obtained within the mean-field calculation is plotted using 
a thick dotted line. The thin dotted line shows the limiting 
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behaviour of x mi for large g, 



lated values of x^ f are shown in Fig. [| together with the 
asymptotic behaviour for large g, x 2 ^ — > 

This behaviour can be explained by thinking in terms 
of an effective potential for the ions V(x). For g < g c , 
the V(x) is a simple harmonic potential. For g > gc , 
it changes into a double-well structure with minima at 
X\ t 2 — ±<?a;mf. In the mean-field approximation no fluc- 
tuations between these minima occur. To go beyond this, 
one needs to include these lattice fluctuations between 
the two minima, which has been considered in Ref. |l2|. 

The magnitude of the lattice fluctuations, (x 2 — (x) 2 ) is 
plotted in Fig. [|. Within the NRG calculation, this quan- 
tity has a clear maximum at a value g* (< g c ). At g c , 
the fluctuations are already significantly reduced. The 
maximum occurs in the crossover region from a single- 
to the double-well potential, where the effective poten- 
tial is broad and shallow. The potential barrier then 
grows rapidly with increasing g. The corresponding fluc- 
tuations in the ME calculation are always small. 

Figure shows the phonon propagator Im d(u) both 
for the Migdal-Eliashberg and the DMFT-NRG calcula- 
tion. The two methods give completely different pictures. 
In the ME approach, the phonon propagator, which for 
9 = consists of a peak at u> = ujq, remains a single peak 
which softens with increasing coupling strength. In con- 
trast to that, the NRG result shows that the main peak 
at lu = ujq broadens, but remains essentially unshifted. In 
addition, a second phonon mode appears at lower energy 
with increasing coupling strength. As g approaches g c 
this mode softens and diverges at g = g c . In the insulat- 
ing phase, only the peak at lu — ujq remains and narrows. 
This behaviour in the insulating phase is to be expected 
as the opening of the gap inhibits any broadening due to 
electron-hole excitations. 

Closer inspection of Figs. |i| and || shows that the soft 
phonon mode develops for g ~ g* , corresponding to the 
maximum in (x 2 — (x) 2 ). This soft mode thus coincides 
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FIG. 5: Phonon propagator o(lu) = — ^Imd(io) for loq = 0.05 
and various coupling strengths g. The upper panel shows the 
Migdal-Eliashberg, the lower panel the NRG results. 



with the build-up of the potential barrier in the effec- 
tive potential of the phonons. From Eq. (||) it follows 
directly that the phonon propagator is closely related 
to the charge susceptibility of the associated impurity 
model. The peak in Im d(u>) has its equivalence in the 
low-energy peak of \c- The existence of such a low-energy 
peak can be expected as follows from the mapping of the 
Holstein onto an attractive Hubbard model. As noted 
before, the physics of the attractive Hubbard model cor- 
respond exactly to those of a repulsive Hubbard model 
with the spin- and charge-channels exchanged. And the 
latter should have a low-energy peak in the spin suscep- 
tibility due to the Kondo physics of its associated impu- 
rity model. The physics of the gap formation, and its 
precursor regime in the Holstein model are thus domi- 
nated by many-body physics. The low-energy feature in 
the phonon propagator has not been predicted before. 

In this letter we have presented the application of the 
dynamical mean-field theory in combination with the nu- 
merical renormalization group to the Holstein model at 
half-filling. This method can be applied to essentially 
all parameter regions of the model. We studied the gap 
formation for small and large phonon frequencies. Gen- 
erally, the gap formation has precursor effects due to 
many-body 'Kondo-like' physics: very strong lattice fluc- 
tuations indicate formation of a double-well potential for 
the ions, and a soft phonon mode emerges due to fluctu- 
ations between the two states of the system. It might be 
possible to observe it experimentally. The framework of 



our method can be extended to contain other local inter- 
actions such as electron-electron interaction of Hubbard- 
type to describe fullerides, or spin-exchange interactions 
as used to describe the manganites. 
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